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Thermodynamic properties of a tetramer ferro-ferro-antiferro-antiferromagnetic Ising-Heisenberg 
bond alternating chain are investigated by the use of an exact mapping transformation technique. 
Exact results for the magnetization, susceptibility and specific heat in the zero as well as nonzero 
magnetic field are presented and discussed in detail. The results obtained from the mapping are 
compared with the relevant experimental data of Cu(3-Clpy)2(N3)2 (3-Clpy=3-Chloropyridine). 
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I. INTRODUCTION 

Quantum behaviour of low-dimensional molecular- 
based magnetic materials has become one of the most 
fascinating topics emerging at the border of condensed 
matter physics, materials science, and inorganic chemis- 
try. In this area, the quantum ferrimagnetic chains 
(QFC) have attracted considerable attention during the 
last few years, because they exhibit a remarkable combi- 
nation of ferromagnetic (F) and antiferromagnetic (AF) 
featuresi^ In an attempt directed toward the synthe- 
sis and design of possible experimental realizations of 
QFC, a magnetostructural analysis of several bimetallic 
assemblies has been accomplished^ since the mixed-spin 
chains afford their most simple and common representa- 
tives. Up to now, the dual aspect of QFC has been ex- 
perimentally confirmed in the [NiCu(pba)(D20)3].2D20 4 
and [MnCu(mal)2 (H20)4]. 2H2O 5 mixed-spin chains. 

In addition to the mixed-spin chains, another class of 
the QFC constitute the bond alternating chains (BAC) 
with an unusual fashion of exchange bonds. From the 
experimental point of view, the bond alternation in the 
ID polymeric assemblies demands at least two struc- 
turally nonequivalent superexchange paths in order to 
get a series of alternating exchange bonds. It should 
be mentioned, however, that bond alternation may also 
arise in a system with an uniform superexchange path- 
way as a result of the spin-Peierls phenomenon, 6 sponta- 
neous dimerization, which occurs when the elastic en- 
ergy increase connected with the lattice distortion is 
lower than the corresponding magnetic energy gain aris- 
ing from the dimerization. Hence, the recent discovery of 



the first inorganic spin-Peierls compounds CuGeOg 7 - and 
a -NaV^Og^ gave rise to a number of theoretical works 
devoted to the spin-1/2 AF-AF BACi2 

Another exciting field in molecular magnetism was 
opened up by Haldane's conjecture^ which has already 
been experimentally verified in several AF nickel-based 
chains (for a review of Haldane gap compounds see Ref. 
fllri . As firstly pointed out by Hida, 12 additional insight 
into the striking properties of the spin-1 AF chain can 
be acquired by analyzing the spin-1/2 F-AF BAC in 
the strong- F-coupling limit 1^ Accordingly, much effort 
has been addressed to prepare polymeric complexes, in 
which the F-AF bond alternation should be realized. At 
present, there exist several copper-based polymeric com- 
pounds that fulfill this requirementjiiiSiiSiiZii&iSiSii (see 
Table I) and Haldane-like behavior has indeed been un- 
doubtedly proved to occur in (IPA)CuCl3. 16 

Of particular interest are also other BAC with more 
peculiar bond alternation, especially with a longer re- 
peating unit of exchange bonds; namely, according to 
the Oshikawa-Yamanaka-Affieck rule*2i one may expect 
the appearance of the magnetization plateau in any sys- 
tem with a longer period of the ground state. Of course, 
this rule represents just a necessary condition for the 
plateau-state formation and does not directly prove its 
existence in any specific model. The theoretical inves- 
tigations focused on the magnetization process of the 
spin-1/2 F-F-AF BAC thus revealed another interesting 
finding - the breakdown of the magnetization plateau. 22 
When the ratio between F and AF coupling constants 
is strong enough in this system, the plateau state disap- 
pears from the magnetization curve. In agreement with 



TABLE I: Several examples of one-dimensional copper-based 
chains with alternating F and AF interactions. 



Chemical formula 


Bond alternation 


Ref. 


[Cu(TIM)]CuCl 4 


F-AF 


114] 


(4-BzPip)CuCl 3 


F-AF 


Il5| 


(IPA)CuCl 3 


F-AF 


I16| 


(DMA)CuCl 3 


F-AF 


[1X1 


[Cu(bipy)(N 3 )2] 


F-AF 


118] 


[Cu2(Me 2 Eten)2(N 3 ) 2 ] (N 3 ) 2 


F-AF 


Il9| 


[Cu 2 (ampy) 2 (N 3 ) 2 ] (N 3 ) 2 


F-AF 


f20l 


3CuCl 2 .2Dx 


F-F-AF 


I23J 


Cu(3-Clpy) 2 (N 3 ) 2 


F-F-AF-AF 


I24J 



Abbreviations: TIM = 2,3,9,10-tetramethyl-l,3,8,10-tetraene- 
cyclo-l,4,8,ll-tetraazatetradecane; 4-BzPip = 4-benzylpipcri- 
dinium (1+) ion; IPA = isopropylammonium (1+) ion; DMA = 
dimcthylammonium (1+) ion; bipy = bipyridine; Dx = 1,4-dioxane; 
Me 2 Eten = N,N-dimethyl-N -ethylethylendiamine; ampy = l-(2- 
aminocthyl)pyrrolidine. 



ledge, there have been reported only a few rigorous re- 
sults for the spin-1/2 F-F-AF-AF BAC obtained by ap- 
plying the exact diagonalization method for a finite-size 
Heisenberg cluster^SSaSLSS and any further more com- 
prehensive studies have not been reported in the litera- 
ture hitherto^ 

The outline of this paper is as follows. In the next 
section, we shall provide the detailed description of the 
model system and then, the basic ideas of transformation 
procedure will be presented. This is followed by the pre- 
sentation of the most interesting results. An exhaustive 
survey of results for several thermodynamic quantities 
in the zero as well as non-zero external field is reviewed 
in Section IIIA, while a comparison with the relevant ex- 
perimental data is included in Section IIIB. Finally, some 
concluding remarks are drawn in Section IV. 



II. MODEL AND METHOD 



this finding, there has not been found any plateau in the 
low-temperature magnetization curve of 3CuCl2-2Dx, 23 
which is regarded as a typical example of the spin-1/2 
F-F-AF BAC with the strong F and weak AF coupling. 

The versatility of the azido ligand in bridging the mag- 
netic ions in various fashions is nicely demonstrated in 
the copper-based compound Cu(3-Clpy)2(N 3 ) 2 r4 here- 
after abbreviated as CCPA. There is a strong ev- 
idence that this compound can be regarded as a 
spin-1/2 tetramer chain with the F-F-AF-AF bond 
alternationi 25 i 26 i 27 ) 28 However, the ID nature of the 
CCPA can be attributed to a sufficient separation be- 
tween chains, which is ensured by the large steric hin- 
drance of the bulky 3-Clpy ligands (see Fig. la). The pe- 
culiar F-F-AF-AF sequence of the exchange bonds arises, 
on the other hand, on account of two kinds of bridging 
fashions of the azido group: the magnetic Cu 2+ ions are 
linked either in a double end-on, or in an end-on and 
end-to-end bridging fashion. It is quite well established 29 
that the end-on bridges are usually associated with the F 
coupling (with exception when the bridging angle is too 
large), while the end-to-end bridges are responsible for 
the AF coupling. The complete magnetic studies of the 
CCPA have been performed by Hagiwara and co-workers: 
for a powder sample high-field magnetization and the 
susceptibility measurements were reported^ while for 
a single-crystal sample the high-field magnetization and, 
susceptibility are known^ as well as the specific hea1i2! 
and electron-spin resonance^ data. 

The primary purpose of this work is to provide a de- 
tailed description to the thermodynamic properties of 
the spin-1/2 F-F-AF-AF BAC by means of a simpli- 
fied Ising-Heisenberg model suggested in our preliminary 
report iSi Due to the simplicity of the proposed Hamilto- 
nian as well as its low-dimensionality, an accurate ana- 
lytical treatment for the complete set of thermodynamic 
quantities can be elaborated within an exact mapping 
transformation technique. 31 To the best of our know- 



A fragment of the CCPA crystal structure is depicted 
in Fig. la. Apparently, there are three nonequivalent po- 
sitions of the Cu 2+ ions due to two different kinds of azido 
bridges. To match the structural situation from the mag- 
netic viewpoint, Fig. lb schematically reproduces the 
magnetic structure of CCPA: the sites interacting purely 
via the F interaction Jp are denoted as the Cul sites, the 
sites coupled via both F as well as AF interaction as the 
Cu2 sites, and finally, the sites coupled purely through 
the AF interaction Jaf are labeled as the Cu3 sites. As 
a consequence of the structural differences, one also has 
to assume various g-factors g\, (72, and 53 at the Cul, 
Cu2, and Cu3 sites, respectively. 

It is worthwhile to say that the supposed magnetic 
structure can also be identified as a ferromagnetic chain, 
the bonds of which are decorated by the AF trimers 
(Fig. lb). Owing to this fact, it is very advisable 
to assume that the F interaction Jp has an Ising-type 
character, whereas the AF interaction Jaf niay have 
the more general form of an anisotropic Heisenberg cou- 
pling (the detailed discussion will be given in Section 
III). Under the circumstances, the model under consid- 
eration can be exactly treated by applying a general- 
ized decoration-iteration mapping transformation^ Ac- 
tually, this mapping procedure has been proved to be 
very useful in investigating several mixed-bond Ising- 
Heisenberg models, some of the present authors already 
obtained within this scheme exact results for the trimer- 
ized Ising-Heisenberg linear chain^ the Ising-Heisenberg 
diamond chain^s, and some decorated Ising-Heisenberg 
planar models. 35 In what follows, we will refer to the 
Cul magnetic sites as to the Ising-type sites, whereas the 
Cu2 and Cu3 sites will be denoted as the Heisenberg-type 
sites. 

Let us write the total Hamiltonian of the spin-1/2 F- 
F-AF-AF Ising-Heisenberg chain comprised of 4A r mag- 
netic sites. By imposing a periodic boundary condition 
(S4JV+1 = Si), the total Hamiltonian of the system takes 



the form: 



;v 



H = JAF 2_^[{Sik-3, S4fc-2)A + (S4fc-2 , S±k-1 ) a] 
fe=l 
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— Jf / „l"4fc-l"4Ai + S ik S ik+1 \ — "1 2-~i 4 



'4k 
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A' 
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4fe-2J 
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fe=l 



where (S.^S^a = A(SfS? + SfSj) + 5?S?, Sf (a = 
x, y, z) marks the spatial components of the spin-1/2 op- 
erator at the zth lattice point and various on-site mag- 
netic fields Bj = gj^sB (j = 1,2,3) have been intro- 
duced in order to distinguish the ^-factors at the Cul, 
Cu2 and Cu3 sites. The first summation accounts for 
the nearest-neighbor AF Heisenberg coupling (Jaf > 0), 
A is the spatial anisotropy in this interaction, and the 
second summation accounts the nearest-neighbor F Ising 
coupling (Jf > 0). Other quantities have the usual 
meaning: \ib is the Bohr magneton and B the external 
magnetic field. 

For convenience, the total Hamiltonian (JJ) can be 
rewritten as a sum of bond Hamiltonians H — ^2 k H k , 
where each bond Hamiltonian Hk involves all the interac- 
tion terms associated with the fcth AF Heisenberg trimcr 
(Fig. lb): 

Hk = JAF[(S4k-3, S4fe-2)A + (S4/C-2, S4fc_i)A] 



MS. 



4k- 



-4^4,k-3 



Sik-iS ik ) 



B ?,Sl k _ 2 



B 2 (Sl k _ x + S z k _ 3 ) - Bx(S z k ^ 4 + S| fc )/2. (2) 



Obviously, the different bond Hamiltonians commute 
with respect to each other, hence, the partition function 
Z can be partially factorized into products of the bond 
partition functions Z k : 



N 



Z = Tr{ 4>8; ... 4J v} Y[ Tr {4fe-3,4fc-2,4fe-l} CXp(~f3H k ) 



Z = Tr 



fe=i 

N 

{4,8, ...4iV} _[_[ %k- 

fe=l 



(3) 



In the above, j3 — l/(fcsT), ks is Boltzmann's 
constant, T is the absolute temperature, the symbol 



Tr 



{4fc-3,4fc-2,4fc-l} 



means the trace over the fcth AF 



Heisenberg trimer, and Tr{4 8j ...an} stands for the trace 
over the spin states of all Ising-type (Cul) spins. To 
proceed further with the calculation, one may introduce 
at the level of the bond partition function Z k an ex- 
tended decoration-iteration transformation by adopting 
the same idea as already discussed in several papers based 
on this mapping scheme (see for instance Refs. |33|-|35|): 



Z k = Tr {4k _ 3Ak _ 2Ak _ 1} eM-PH k ) = exp[/3B 1 (5| fc _ 4 + Sf fe )/2j <^ exp[/3(J j4F + J+)/6] ^ exp(/3a; n ) + (4) 

^ n=l 

exp[(3(J AF - J+)/6] J2 e MPVn) + 2exp(- / 3J j4F /2) cosh(/3J+/2) j = Aexp[f3RS z 4k _ 4 S z 4k + f3H (S z k _ 4 + S z k )/2}. 
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Since a complete explicit rendering of the decoration- equations that are given as follows: 
iteration mapping is in this case rather intricate, we have 
defined in the transformation Eq. @ several functions 
in order to write it in a more appropriate form: 



Jf 


= Jp(S 4 z k _ 4 + S% k ) + 2B 2 + B 3 , 


(5) 


Jf 


= J F (S z k _ i + Sl k ) + 2(B 2 - B 3 ), 


(6) 


JT 


= JF\S Ak _ i — S ik ), 


(7) 



and the terms x n and y n denote the roots of two cubic 



x n = ±2i/Pl cos[$i + (n - 1)2tt/3], (n = 1, 2, 3), (8) 
Pi = [(Jaf + Jf)/&f + [(Jf) 2 + 2(J ylF A) 2 ]/12, (9) 
[(Jaf + Jf)/&} 3 + 

+ [(Jaf + J 2 + )/6][(J AF A) 2 - (Jf) 2 ]/4, (10) 

' ' - (ID 



Q 



$i = - arctan (JP? - Q\ /Qj. 



and 



y n = ±2v^cos[$ 2 + (n - 1)2tt/3], (n = 1, 2, 3), (12) 
P 2 = [(Jaf - J 2 + )/6] 2 + [(./r) 2 + 2(J j4F A) 2 ]/12, (13) 
Q 2 = [(Jaf - J 2 + )/6] 3 + 

+ [(J AF - J+)/6][(J AF A) 2 - (Jf) 2 ]/4, (14) 



*2 = r arctanf \AR? - Q 2 



i/Qi 



(15) 



Note that the signs of x n and j/„ are unambiguously de- 
termined by the signs of the expressions Q% and Q 2 , re- 
spectively. 

It should be emphasized that the mapping parameters 
A, R and Hq are " self-consistently" given by the trans- 
formation Eq. Q, which must be valid for any combi- 
nation of spin states of the Sik-i and S 1 ^ Ising spins. 
In consequence of that, the mapping parameters can be 
obtained following the standard procedure 31 - 33 ' 34,35 from 
the expressions: 

A 4 = V 1 V 2 Vi, pR = \n(V l V 2 /V£), /3H =]n(y 1 /V 2 ), 

(16) 

where the functions V% , V 2 and V 3 have a physical mean- 
ing of the bond partition function Q, when taking into 
account all the particular spin combinations of the S^k-i 
and Sik Ising spins. These, however, modify merely the 
effective coupling constants ©-0 implicitly contained 
in Eq. ||3J). Thus, the functions V\, V 2 , and V 3 are defi- 
nitely determined by this set of expressions: 

V 1 = Z k /±J+ = J F + 2B 2 + B 3l Jf = 0, 

J+ = J F + 2(B 2 - B 3 ); (17) 

V 2 = Z k , if: J+ = -J F + 2B 2 + B 3 , Ji = 0, 

4 = -J F + 2{B 2 -B 3 )- (18) 

V 3 = Z k , if: J+ = 2B 2 + B 3 , Ji = J F , 



J+ = 2{B 2 - B 3 



(19) 



When substituting Eq. (gj into Eq. 0, the trans- 
formation Eq. JIJ maps the original F-F-AF-AF Ising- 
Heisenberg BAC on the uniform spin-1/2 Ising chain with 
the effective nearest-neighbor exchange coupling R and 
the magnetic field Hq. As a result, the partition function 
Z of the Ising-Heisenberg BAC can be directly related 
to the partition function Zq of the corresponding Ising 
chain: 



Z = A N Z (P,R,H ). 



(20) 



Certainly, similar mapping relations can be established 
also for other thermodynamic quantities. For instance, 
the Gibbs free energy G of the Ising-Heisenberg BAC can 
be calculated from the relevant expression of the Gibbs 
free energy Go of the corresponding spin-1/2 Ising chain: 



G = G (/3,R,H ) -Nk B Tln(A). 



(21) 



Since an exact solution for the spin-1/2 Ising chain was 
known a long time ago^ the above equation can serve as 



the basic generating equation from which all thermody- 
namic quantities can be extracted. Here, we shall restrict 
ourselves to the analysis of the magnetization, suscepti- 
bility and specific heat. For the spin only value of the 
on-site magnetization at the magnetically nonequivalent 
Cul, Cu2 and Cu3 positions, we shall introduce this sim- 
ple notation: 

m 1 = {Si k ), ma = <&_!>, m 3 = {Sl k _ 2 ), (22) 

where the symbol (...) stands for a standard canonical 
average over the ensemble defined by the Hamiltonian 
(Q. In view of this notation, the total magnetization 
normalized per one Cu 2+ ion can be expressed as follows: 
to = H B {gimi + 2g 2 m 2 + g 3 m 3 )/4. 

Finally, we briefly mention the basic thermodynamic 
relations from which all the analyzed quantities have 
been calculated after straightforward but a little bit 
lcngthly calculations. The on-site magnetization can be 
obtained by differentiating the Gibbs free energy with 
respect to the particular magnetic fields: 



TOl 



1 dG 
NdB^ 



,m 2 



1 dG 

2N8B- 



■,m 3 



1 dG 
NOBl 



(23) 



while the susceptibility and specific heat have been ob- 
tained as the second derivatives of the Gibbs potential 
using the standard thermodynamic relations: 



X = 



d 2 G 
dB 2 ' 



and 



C = -T 



d 2 G 
dT 2 ' 



(24) 



It should be stressed that the explicit form for these quan- 
tities is too cumbersome to write it down here; however, 
it can be obtained from the authors on request. 



III. RESULTS AND DISCUSSION 

Before proceeding to a discussion of the most inter- 
esting results, the model reliability should be checked in 
connection with its possible application to interpret the 
experimental data on CCPA, because a danger of over- 
interpretation is inherent in any approximation. At first 
glance, we have made in our model a very crude conjec- 
ture in that the F interaction J F was approximated by 
an Ising-type coupling although all Cu 2+ ions are nearly 
isotropic, and whence, the Heisenberg interaction would 
be more appropriate. With regard to this, Fig. 2 il- 
lustrates a comparison between the on-site magnetiza- 
tion of the Ising-Heisenberg BAC and that of the cor- 
responding pure Heisenberg BAC obtained from the ex- 
act diagonalization method for a finite-size cluster of 12 
sitesiSi It should be pointed out that the low-temperature 
{ksT/ Jaf — 0.05) magnetization curve provides the best 
independent test of the model reliability, since it reflects 
magnetization near the ground state, where the clear- 
est manifestation of the quantum fluctuations should be 
expected to occur. 



It is quite surprising that there is such an excellent 
agreement between the two theories: the total magneti- 
zation (for clarity not shown here, but see for instance 
Fig. 4) exhibits a steep increase from zero field followed 
by a magnetization plateau and finally, there appears a 
second steep increase near the transition field toward a 
fully polarized state. Strictly speaking, there is no real 
phase transition at any finite temperature; however, the 
low-temperature magnetization curve shed light on what 
happens in the ground state: the plateau state should 
reflect the ground-state phase and the field-induced tran- 
sition to the fully saturated phase takes indeed place in 
the zero-temperature limit. Moreover, the magnetic or- 
der at the plateau state has a typical feature of quantum 
fcrrimagnet, in fact, one finds here a substantial quan- 
tum reduction of the magnetization m^ and m^ at the 
Cu2 and Cu3 sites even though the magnetization m\ 
(Cul site) retains its saturation value. It means, among 
other matters, that the spin deviations cannot propagate 
through these sites and hence, they are strictly localized 
within the AF trimers. This explains also the remarkable 
agreement between both the theories, since the quantum 
fluctuations are in our model artificially restricted to the 
AF trimers because of the presence of the Ising spins. 

The most obvious difference between the magnetiza- 
tion curve of the Ising-Heisenberg BAC and the pure 
Heisenberg BAC can be thus found in the vicinity of 
the zero field, where the on-site magnetization of the 
former reaches their plateau values more rapidly. Ap- 
parently, this distinction can be attributed to the more 
"susceptible" character of the Ising spins, whereby their 
Heisenberg counterparts achieve the plateau values more 
rapidly on behalf of the local field produced just by the 
Ising spins. With this in mind, small systematic dis- 
crepancies should be expected to occur especially in the 
region where both B — » and T — » 0, but, without loss 
of the qualitative agreement. 

Why the Cul sites represent a "barrier" for the spin 
deviations also in the spin-1/2 Heisenberg F-F-AF-AF 
BAC still remains an open question. Although the origin 
of this outstanding feature could be naively understood 
as a general feature of the spins coupled by the F interac- 
tions only, this suggestion is in an obvious contradiction 
with Hida's results for the spin-1/2 F-F-AF BAC^ As 
a matter of fact, the total magnetization of this system 
varies smoothly with increasing magnetic field even in the 
ground state, which implies that the spin deviations are 
delocalized over the whole chain and thus, there is a cer- 
tain quantum reduction also at purely ferromagnetically 
coupled sites. 



A. Survey of theoretical results 

In this part, an extensive survey of theoretical results 
will be presented in order to enable deep insight into 
the magnetic behavior of the system under investigation. 
For the sake of simplicity, we shall firstly suppose equal 



g-factors at the Cul, Cu2, and Cu3 positions, i. e. g\ — 

92 = 93 = 9- 

We start our discussion with the analysis of the ground 
state. Below some critical external magnetic field, the 
ground state exhibits an interesting fcrrimagnetic order, 
which can be characterized by the following values of the 
ground-state magnetization: 



1 



1 



mi 



m 3 



2' m2 = 4 



1 + Jf/Jaf 



■ 70" 

Jf/Jaf 



2 X /(1 + J F /J^)2+8A 



Jf/Jaf) 2 + 8A* 
m 1 
2 



^- = 5. (25) 



where m s labels the saturation magnetization normalized 
per one Cu 2+ ion. As already mentioned in the preceding 
part, the magnetization m\ retains in the ground state 
its saturation value in contrast to the magnetization m-i 
and 7713, which show obvious quantum reduction. It is 
quite evident that the stronger the exchange anisotropy 
A, the greater the reduction of the magnetization to 2 
and 7713. On the other hand, the quantum reduction of 
the magnetization completely vanishes in the Ising limit 
A — » 0, when the ground state exhibits a "classical" ferri- 
magnetic order TTIT ■ ■ ■ (CulCu2Cu3Cu2...) of a rather 
trivial nature. 

As expected, the ferrimagnetic system undergoes a 
field-induced metamagnetic transition toward the fully 
saturated phase under the certain magnetic field. It can 
be readily proved that the transition field B t is given 
under the condition: 



gfj, B B t /JAF = -t[3- Jf/Jaf 



V(l + Jf I Jaf) 2 + 8A2 . (26) 



For illustration, Fig. 3 displays the ground-state phase 
diagram in the Jp — B plane for several values of the ex- 
change anisotropy A. As one can see, increasing strength 
of the F coupling Jp generally reinforces the gradual 
decline of the transition field. Thereby, the pure Ising 
limit (A = 0) represents the only exceptional case when 
the magnitude of the transition field does not change as 
the ratio Jf/Jaf varies. It can be easily understood 
that the metamagnetic transition arises in this particu- 
lar case from a single flip of the central (Cu3) spin of 
each AF trimer, which occurs when the exchange energy 
of Jaf is thoroughly balanced by the external magnetic 
field (gnBB t = Jaf)- 

The situation becomes much more complex on consid- 
ering the nonzero exchange anisotropy A. Due to the 
quantum fluctuations, the spin oriented in an opposite 
direction with respect to the external field is no longer 
rigidly connected with the central spin of the AF trimer, 
but is collectively held by the entire trimer. In other 
words, the reversed spin is delocalized over the whole AF 
trimer and consequently there is a non-zero probability 
that the side (Cu2) spins of the trimer are aligned op- 
posite to the field direction. Naturally, this must lead 



to the enhancement of the exchange energy between the 
ferromagnetically coupled Cu2 side spins and the fully 
polarized Cul spins. In consequence of that, it is quite 
conspicuous that the suppression of the transition field 
can be explained in terms of the energetic destabilization 
of the quantum ferrimagnetic order, which occurs when 
Jf/Jaf increases. 

Now, let us turn our attention to the magnetic behav- 
ior at finite temperatures. In Fig. 4 we plot the total 
magnetization scaled in gfis units against the dimen- 
sionless magnetic field (gfisB/ ' Jaf) for a few temper- 
atures. It should be pointed out that the magnetization 
curve starts at any finite temperature from zero accord- 
ing to the one-dimensional character of the spin system. 
The sharp stepwise magnetization curve observable in the 
zero-temperature limit is, however, gradually smeared 
out as the temperature is raised from zero. Actually, 
the conversion toward the fully saturated phase does not 
occur merely at one precise value of the transition field, 
but is smudged over a finite range of the fields. In addi- 
tion, the temperature-induced fluctuations greatly shrink 
also the width of the magnetization plateau and above a 
certain temperature, the plateau state completely disap- 
pears from the magnetization curve. By any means, Fig. 
4 provides further convincing evidence that the observed 
magnetization plateau emerges, under the assumption of 
the uniform ^-factors, exactly at one-half of the satura- 
tion magnetization [see also Eq. 1)25(1], 

The variation of the field-induced magnetization with 
the temperature (Fig. 5) is also interesting on account 
of the magnetically ordered ground state. At low mag- 
netic fields, the total magnetization tends abruptly to 
zero as the temperature increases. Nevertheless, one 
finds here also a noticeable temperature- induced increase 
of the magnetization within a range of moderate fields, 
when the total magnetization shows a broad maximum 
resulting from a vigorous thermal excitations of the Cu3 
spins and smaller thermal excitations of the Cu2 spins. 
However, when the magnetic field exceeds the transition 
field given by the condition IJ26II , the total magnetization 
monotonically decreases with increasing temperature as 
it already starts from its maximum value. 

The thermal dependence of the zero-field susceptibility 
times temperature (xT) data is displayed in Fig. 6. In- 
terestingly, the \T product exhibits a round minimum 
upon cooling and then diverges under further tempe- 
rature suppression. As already discussed in several pa- 
pers concerned with the QFC^ a thermal dependence 
of this type reveals a ferromagnetic-to-antiferromagnetic 
crossover. Really, the temperature variation of the yT 
data is in general a monotonically decreasing function for 
ferromagnets, but monotonically increasing for antifcrro- 
magnets. In this respect, the marked low-temperature di- 
vergence of the susceptibility emerges markedly because 
of the gapless F excitations from the magnetically or- 
dered ground state. The position of the round minimum, 
on the other hand, designates the temperature above 
which the gapped excitations of the AF nature overhelm 



the gapless excitations originating from Jp. In accor- 
dance with the above statement, the round minimum 
flattens as the ratio Jf/Jaf strengthens and simultane- 
ously its position is shifted toward higher temperatures. 

As far as the magnetic susceptibility at nonzero fields 
is concerned (Fig. 7), the \F always initiates from zero 
because of the energy gap opened by the magnetic field. 
However, the thermal dependence of the \F data shows 
at weak fields a dramatic increase until it reaches a sharp 
maximum, which is followed by the familiar round min- 
imum of the same origin, as already discussed by the 
zero-field susceptibility. The appearance of the addi- 
tional sharp peak can obviously be related to a ther- 
mal instability of the magnetically ordered system, be- 
cause already a small temperature change necessitates a 
huge variation of the magnetization at low fields (see Fig. 
5). As a matter of fact, the yT product rises steadily 
with the temperature at sufficiently strong fields, since 
the thermal fluctuations are not strong enough to in- 
duce a simultaneous excitation of a large number of spins. 
Thus, an interesting thermal dependence of the suscep- 
tibility appears for strong magnetic fields just around 
the transition to the fully polarized state. Under these 
circumstances, the T-\T plot exhibits a rapid increase 
over a small temperature range, which is subsequently 
followed by a narrow plateau that continuously passes 
into a slowly repeating increase of the susceptibility (see 
the curve g^isB / Jaf — 1-4). 

To receive a complete picture of the thermodynam- 
ics, the thermal variations of the specific heat are plot- 
ted in Figs. 8 and 9. When considering the zero-field 
case (Fig. 8), the specific heat displays as a function 
of the temperature a remarkable double-peak structure. 
There are strong indications that the first sharper peak 
observed at lower temperature originates from the F ex- 
citations: the peak becomes wider as the ratio Jf/Jaf 
increases and is shifted to higher temperatures notwith- 
standing its almost constant height. Under the assump- 
tion of dominant F coupling constant {Jf/Jaf > 1-0), 
the low-temperature peak is to a large extent superim- 
posed on the second broader maximum that occurs at a 
little bit higher temperature. Obviously, the round high- 
temperature maximum can be thought of as the usual 
Schottky-type maximum, which has a tendency to be en- 
hanced in magnitude with increasing Jf/Jaf- When 
comparing the results displayed in Fig. 8 with those 
of the pure Heisenberg F-F-AF-AF BAC obtained us- 
ing the exact diagonalization method^ or the transfer- 
matrix renormalization group 3 ^, an excellent agreement 
is found as far as the high-temperature maxima of both 
these models are concerned. The only difference thus 
rests in the height of the low-temperature peak which 
is, due to the Ising approximation of Jp, roughly three- 
times higher for the Ising-Heisenberg BAC compared to 
that of the pure Heisenberg BAC (see for instance the 
insert of Fig. 4 depicted in Ref. I32T) . 

The situation becomes even much more interesting on 
applying a magnetic field (Fig. 9). The effect of a 



small magnetic field is to increase the height of the low- 
temperature peak and to move it toward higher temper- 
atures (see the case g/j, B B/JAF — 0.1). This result is 
taken to mean that under a certain field both maxima 
coalesce and consequently, the specific heat exhibits a 
single nonroundcd maximum as shown for gfJ, B B / J af — 
0.2. By an additional increase of the field strength, the 
specific heat curve gradually loses its irregular profile 
and the overall trend is that the Schottky-type maxi- 
mum drops in magnitude and moves to higher tempera- 
tures. Apart from this rather trivial finding, the double- 
peak specific heat curve can be recovered for the fields 
close to the transition to the fully saturated phase (e.g. 
gii B Bj Jaf = 1-4). It should be mentioned, however, 
that the height of the low-temperature peak is in this 
case considerably smaller than that of the zero-field spe- 
cific heat. Although the external field spreads also this 



J 



low-temperature peak until it completely merges with the 
Schottky-type maximum, the change of the peak posi- 
tion occurs strikingly without any significant change of 
the peak height (for clarity, this effect is not shown here). 
Finally, it should be stressed that the observed behavior 
of the specific heat is not a generic feature of the spe- 
cial class of the Ising-Heisenbcrg BAC, nevertheless, it 
has already been recognized also in the pure Ising BAC 
(compare with Figs. 4 and 8 from Ref. l37l) . 

We conclude our survey of the thermodynamic pro- 
perties by considering the nonuniformity of the g- 
factors at the structurally and therefore also magneti- 
cally nonequivalent Cul, Cu2 and Cu3 positions (i.e. 
,<?i 7^ .92 ^ g 3 )- Let us firstly focus on the ground-state 
behavior. The most interesting finding to emerge here is 
that the ground-state magnetization depend on the field 
strength, whenever g 2 7^ 33^ 



mi = 



1 
m 3 = --- 



1 , 
mo = -< 1 



1 + Jf/Jaf + 2fi B B(g 2 - g-i)/JAF 
4 l~ ' VI 1 + Jf/Jaf + 2 l i B B{g 2 - g 3 )/J AF } 2 + 8A* 
1 + J F J J A F + 2fi B B{g 2 - g 3 )/J A F 



2 y/[\ + J F /J AF + 2n B B(g 2 - g 3 )/J AF ] 2 + 8A 2 ' 

m fli +ff2 ff2 -ff3 1 + Jf/Jaf + 2fi B B(g 2 - g 3 )/J AF 

m s g\ + 2g 2 + g 3 g x + 2g 2 + g 3 ^[1 + J F /J AF + 2^i B B{g 2 - g 3 )/J AF } 2 + 8A 2 



(27) 



r 



It is quite apparent from Eqs. I|27[) that the field de- 
pendence of the total magnetization comes from the cor- 
responding field variations of the on-site magnetization 
m 2 and m 3 . Another noticeable feature to observe here 
is that the magnetic field suppresses (raises) the quan- 
tum reduction of the magnetization m 2 and m 3 as far as 
.92 > 33 (.92 < 33)- Since the condition 2\i B B{g 2 - g 3 ) < 
Jf + Jaf holds for most of the experimentally accessi- 
ble fields and 3-factors, the total magnetization should 
be nearly linearly dependent on the magnetic field with 
a linear term proportional to the difference Sg = g 2 — g 3 . 
Naturally, the greater the magnitude of the external field 
and 8g, the stronger deviations from the linearity should 
emerge. 

Some typical examples of the ground-state magnetiza- 
tion curves are depicted in Fig. 10 under the assumption 
.92 > .91 = 33 (Fig- 10a) and 33 > 31 = g 2 (Fig. 10b). 
These results can serve as evidence that the magneti- 
zation curve does not show an exact plateau, but rises 
steadily with the external field. It is also worthy to note 
that these "magnetization plateaus with a finite slope" 
(when plotted in the full range with respect to the satu- 
ration magnetization, they are hardly discernible from 
the exact plateaus within the reasonable values of the 3- 
factors) do not occur precisely at half of the saturation 
magnetization {rn/m s — 1/2), but they are shifted to 
higher values if g 2 > g\ — 33 (Fig. 10a) and respectively, 



to the smaller values if 33 > g\ — g 2 (Fig. 10b). 

To enable an independent check of the field dependence 
of the total magnetization, the temperature variation of 
the susceptibility has been examined in detail. The typi- 
cal thermal dependences of the susceptibility are plotted 
in Fig. 11 for g 3 > gi — g 2 and some selected exter- 
nal fields. Generally, the susceptibility versus temper- 
ature plot can be characterized by a round maximum, 
which flattens and shifts to higher temperatures when 
increasing the external field. The low-temperature part 
of the susceptibility is displayed on an enlarged scale in 
the inset. It turns out that the susceptibility does not 
completely vanish as the temperature approaches zero; 
nevertheless, it tends toward a small but finite value. 
This result provides a confirmation of the striking field- 
dependent magnetization and moreover, it also clearly 
demonstrates a gapless excitation spectrum above the 
ground state. 



B. Comparison with the experimental data 

At this stage, the results obtained from the exact map- 
ping on the Ising-Heisenberg BAC will be compared with 
the available experimental data for a single-crystal sam- 
ple of the CCPA. It should be mentioned here that all the 
experimental data used in our further analysis have al- 
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ready been reported in earlier publications 26,27 to which 
the interested reader is referred for more experimental 
details. 

In order to fit the experimental data of the CCPA, 
we have to consider the six fitting parameters: the ex- 
change constant Jaf/^b, the exchange anisotropy A, the 
alternation ratio Jf/Jaf and the three g-factors g\ 1 gi 
and 53. Note that the best fit of the magnetization and 
the susceptibility data obtained under the restriction of 
equal ^-factors {g\ = 52 = 53) as well as isotropic ex- 
change Jaf (A = 1.0), have already been published by 
the present authors in our preliminary report (see Fig. 3 
in Ref. l30f) . Unless specifically mentioned, the exchange 
anisotropy will again be fixed to the purely isotropic case 
A = 1.0, because the experimental measurements on the 
single-crystal sample do not reveal any significant spatial 
dependence: the spatial directions parallel and perpen- 
dicular to the chain axis seem to be almost equivalent. 26 

In Fig. 12, the low-temperature magnetization curve of 
CCPA measured in the pulsed (Fig. 12a) and static (Fig. 
12b) magnetic fields is compared with the relevant the- 
oretical prediction (the fitting parameters are indicated 
in the figures). It turns out, however, that the average 
g-factor must be approximately g av = {g\ + 2^2 +33) /4 ~ 
2.15 in order to get the correct saturation magnetization 
(Fig. 12a). Even under this confinement, several com- 
binations of the g-factors provide almost the same fit to 
the experimental magnetization curve. For illustrative 
purposes, we choose the one with the smallest value of 
the gi-factor to lower the susceptibility of the Cul spins, 
because the most obvious discrepancies between the the- 
ory and the experiment emerge near zero field due to 
the more "susceptible" character of the Ising spins, as 
already discussed earlier. 

Unfortunately, it is beyond experimental verification 
to find out from the existing data set, whether the mag- 
netization increase observed within the plateau region 
appears exclusively on account of the finite-temperature 
effect, or is partially caused by the inequality 32 7^ .93- 
Although the former situation cannot be definitely ruled 
out, it is rather more acceptable that there is at least 
a small contribution to the magnetization increase also 
from the difference between the ^-factors. Nonetheless, 
the final conclusion is complicated by different factors: 
the magnetization curves in the static fields (Fig. 12b) 
are measured at higher temperatures and become quite 
noisy above 15 T, while the available magnetization curve 
in the pulsed fields (Fig. 12a) show a relatively large hys- 
teresis probably caused by the magnetocaloric effect. 26 

The thermal dependence of the zero-field susceptibility 
times temperature data is illustrated in Fig. 13. To ob- 
tain a reasonable accordance with the experimental mea- 
surement, slightly higher g- factors must be taken into 
consideration in comparison with the ones acquired by 
the fitting of the magnetization. Despite this fact, a 
quite impressive accordance can be achieved by rescaling 
the average g-factor at about 4-5% of its magnetization 
value g av ^2.15. As expected, the most apparent dis- 



agreement appears in the descending tail of the yT curve 
nearby the zero temperature, where the susceptibility di- 
vergence is reminiscent of the ferromagnetic Ising-type 
divergence rather than the Heisenberg-type divergence. 

Finally, Fig. 14 shows the magnetic part of the specific 
heat as a function of the temperature for several magnetic 
fields. The corresponding theoretical predictions are also 
included; nevertheless, the specific heat can merely be fit- 
ted with a rather unacceptable precision even though the 
results agree at least qualitatively. In order to situate the 
striking low-temperature peak around 0.5 K, however, a 
drastic reduction of the Jf and Jaf exchange parameters 
was carried out. Even under this condition, the predicted 
peak is much more robust than observed experimentally 
and this feature remains valid also for any other combi- 
nation of the Jf and Jaf coupling constants. Actually, 
the height of the low-temperature peak does not signifi- 
cantly vary with the ratio Jf/Jaf as already shown in 
Fig. 8. 

It is quite obvious that the inconsistency in the fitting 
set for the specific heat and, respectively, the magnetiza- 
tion and susceptibility, indicates some insufficiency in the 
considered model based on the simple concept of the F-F- 
AF-AF BAC. Notice that the basic problem that prevents 
obtaining an unique fitting set for all these quantities con- 
sists in reproducing the striking low-temperature peak of 
the specific heat, which cannot be fitted within a reason- 
able accord neither under the assumption of the pure 
Heisenberg F-F-AF-AF BACiSL 3 ^ A possible explana- 
tion for this peculiar inconsistency has recently been sug- 
gested by Lu et aL— , who found a strong indication that 
the striking low-temperature peak comes from the con- 
tribution of magnetic impurities. Until a higher-quality 
sample of CCPA is prepared and remeasured, however, 
one cannot exclude the other possibility that this anoma- 
lous peak originates from the neglected higher-order in- 
teraction terms in the proposed Hamiltonian. Our de- 
tailed analysis shows that neither the nonuniformity of 
g-i actors nor the exchange anisotropy of the Jaf (A) cou- 
pling constant, can resolve this inconsistency. 



IV. CONCLUDING REMARKS 

In the present article, the magnetic properties of the 
spin-1/2 Ising-Heisenberg chain with regular F-F-AF-AF 
bond alternation have been investigated within the ex- 
act mapping transformation method. The results ob- 
tained from the mapping procedure are compared with 
the results of the corresponding Heisenberg BAC 26 and 
with the experimental data of the CCPA compound^SiS 7 . 
which is regarded to be an experimental realization of 
the spin-1/2 F-F-AF-AF BAC. 

It should be emphasized that all characteristic quan- 
tum features observed in the spin-1/2 Heisenberg F-F- 
AF-AF BAG— £i£&2£ are still present also in our simpli- 
fied Ising-Heisenberg version: one finds here a substantial 
quantum reduction of the ground-state magnetization mi 



and rri3 at the Cu2 and Cu3 sites in addition to the fully 
saturated ground-state magnetization mi at the Cul site. 
Perhaps, the most striking finding stemming from our 
study is that there appears a rather unusual change from 
the gapful to gapless excitation spectrum when convert- 
ing the uniform g-factors at structurally non-equivalent 
positions to nonuniform ones. With regard to this, the 
magnetization curve exhibits a real plateau under the 
condition (72 = 53 only; for any other case the magneti- 
zation within the quantum ferrimagnetic phase (plateau 
state) shows a weak field dependence. 

The success of the simplified Ising-Heisenberg model in 
reproducing the measured data for CCPA is also quite re- 
markable. With the exception of the specific heat, where 
the results agree at the qualitative level only, other ther- 
modynamic quantities are relatively well reproduced. In 
fact, small systematic errors occur merely in the low- 
temperature and weak-field region, where the approxima- 
tion of the ferromagnetic coupling Jp by the Ising-type 
interaction plays the most essential role. It is worthwhile 
to mention, moreover, that the discrepancy found in the 
specific heat fit cannot be resolved even if a pure Heisen- 
berg F-F-AF-AF BAC is considered. Actually, very re- 
cent work of Lu et al£2- provided a strong indication that 
the striking low-temperature peak of specific heat is be- 
ing of extrinsic origin probably caused by the presence 
of magnetic impurities. However, further experimental 
measurements with high-quality samples of CCPA are 
required to clarify this issue. Finally, one should notice 



that the simple concept of the Ising-Heisenberg F-F-AF- 
AF BAC presented here can be rather straightforwardly 
extended to account for second-neighbor interactions, the 
antisymetric Dzyaloshinsky-Moriya interaction, and/or 
multispin interactions, which cannot be ruled out as a 
potential cause of the observed discrepancies until the 
influence of magnetic impurities will experimentally be 
confirmed. 
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Figure captions 

Fig. 1 The fragment of the CCPA structure. Fig. la 
shows a full details with the azido (N3) bridges as 
well as the bulky 3-Clpy ligands, Fig. lb schemati- 
cally reproduces the considered magnetic structure 
only. There are three non-equivalent positions of 
the Cu 2+ ions: the Cul sites coupled purely by the 
F interaction (black circles), the Cu2 sites coupled 
via both F and AF interactions (grey circles) and 
lastly, the Cu3 sites coupled purely by the AF inter- 
action (light-grey circles). The rectangles designate 
the AF trimers. 

Fig. 2 The comparison between the on-site magnetization 
of the Ising-Heisenberg BAC and the corresponding 
pure Heisenberg BAC obtained by the exact diag- 
onalization method as described in Ref. 26. For 
details see the text. 

Fig. 3 The ground-state phase diagram in the Jp-B plane 
obtained for several exchange anisotropies A under 
the assumption of the uniform ^-factors, i. e. g\ = 
92 = .93 = 9- 

Fig. 4 The total magnetization scaled in the g/is units 
as a function of the dimensionless magnetic field 
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gfisB/JAF for Jf/Jaf — 0.5, A = 1.0 and several 
dimcnsionless temperatures ksT ' j ' J af- 

Fig. 5 Some typical temperature variation of the total 
magnetization obtained for Jf/Jaf — 0.5, A = 1.0 
and several magnetic fields. 

Fig. 6 The thermal dependence of the zero-field sus- 
ceptibility times temperature data scaled in the 
NA(g/is) 2 /fcs units (Na - Avogadro's number) for 
A = 1.0 and some typical values of the Jf/Jaf- 

Fig. 7 The temperature dependence of the \F product as 
a function of the field strength when Jf/Jaf = 0.5 
and A = 1.0. 

Fig. 8 The zero- field specific heat as a function of the tem- 
perature for some typical values of the Jf/Jaf and 
A = 1.0. The specific heat is scaled in multipliers of 
the universal gas constant R = 8.314 J.K _1 .mol _1 . 

Fig. 9 The thermal dependence of the specific heat by var- 
ious magnitudes of the external magnetic field when 
Jf/Jaf = 0.5 and A = 1.0. 

Fig. 10 The field variation of the ground-state magnetiza- 
tion normalized with respect to its saturation value 
by assuming the non-uniformity of the ^-factors. 



Fig. 10a (Fig. 10b) shows the situation when 
32 > 3i =53 (33 > 3i =32)- 

Fig. 11 The temperature dependence of the susceptibil- 
ity when varying the external field strength, and 
Jf/Jaf = 0.5, A = 1.0, g x = g 2 = 2.0 and 
93 = 2.25. The insert displays the low-temperature 
susceptibility in enlargened scala. 

Fig. 12 The low-temperature magnetization curves mea- 
sured in the pulsed (Fig. 12a) and static (Fig. 12b) 
magnetic fields together with the corresponding 
theoretical prediction obtained for this fitting set 
of parameters: Jaf/^b = 28.2 K, Jf/Jaf = 0.4, 
A = 1.0, g x = 2.10, 52 = 2.14 and g 3 = 2.19. The 
hysteresis observed in the experimental curves near 
the zero and saturation fields is probably caused by 
the magnetocaloric effect. 

Fig. 13 The zero-field susceptibility times temperature 
data and the corresponding theoretical prediction 
(the fitting parameters arc indicated in the figure). 

Fig. 14 The thermal dependence of the specific heat mea- 
sured at B — 0.0, 0.5 and 1.0 T. The corresponding 
theoretical predictions are also included. 
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